Constrained Brownian motion as thermodiffusion 
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We propose a Langevin equation for systems in an environment with nonuniform temperature. 
At odds with an older proposal, ours admits a locally Maxwellian steady state, local equipartition 
holds and for detailed-balanced (reversible) systems statistical and physical entropies coincide. We 
describe its thermodynamics, which entails a generalized version of the First Law and Clausius's 
characterization of reversibility. Finally, we show that a Brownian particle constrained into a smooth 
curve behaves according to our equation, as if experiencing nonuniform temperature. 
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It is well known that noise drives particles along the 
temperature gradient, giving rise to thermophoresis, or 
thermodiffusion. "Pebbles in a driveway accumulate on 
the side" wrote Landauer suggestively "In the driveway 
they are agitated (hot region). They are left undisturbed 
on the side (cold region)" [3. He had long sharpened his 
intuition on the problem [2|, and came to argue that, in 
the light of the equipartition theorem for the root mean 
square velocity v ~ VT, detailed balancing of the cur- 
rents PhVh and p c v c , respectively out of a hot region and 
of a cold one, implies a spatial density p ~ 1/y/T. He 
admitted though that this all-too-simple argument did 
not do justice to more complicated situations where the 
geometry of the mean free path changes with temper- 
ature, and that other proposals were valuable, in par- 
ticular Van Kampen's p ~ 1/T, which followed from the 
analysis of a Langevin equation with state-dependent dif- 
fusion coefficient [3|. Since then, Van Kampen's theory 
has been used by various authors to discuss the stochas- 
tic energetics of nonuniform temperature systems [ij (see 
also Sekimoto's monograph [fj|, §1.3.2.1, §4.1.2.2]), to de- 
rive the Soret coefficient Q, and more recently to find 
an anomaly between overdamped thermodynamics and 
thermodynamics with finite, but small, inertia [7J. 

However, when applied to reversible (i.e. detailed bal- 
anced) systems, Van Kampen's theory displays some odd 
peculiarities. The steady state is not locally Maxwellian 
(Local Thermal Equilibrium). Equipartition fails: Tem- 
perature is not the average kinetic energy. The conven- 
tional definition of heat does not allow to match statisti- 
cal and physical notions of entropy, a cornerstone princi- 
ple of Equilibrium Statistical Mechanics. Moreover, gen- 
eralization to rt-dimensional systems with non-diagonal 
diffusion matrix leads to nonequilibrium steady states 
even in the absence of external forces Q. All these ef- 
fects, if a bit controversial from a foundational perspec- 
tive, might well occur and are of great interest: After all 
equipartition and Gaussianity should be ascribed to ideal 
situations. However, so far we lack a basic ideal model. 



1. Scope and overview 

In this work we enforce the principle of Local Ther- 
mal Equilibrium to obtain a new and different theory 
that marries well with the thermodynamics of nonuni- 
form temperature systems and panders to Landauer's 
intuitions. We first derive the relevant equations and 
discuss their dynamical properties and time scales. We 
then define proper equilibrium thermodynamic poten- 
tials, discuss conditions for the occurrence of equilib- 
rium and dwell into thermodiffusion in nonequilibrium 
regimes. We already pinpoint that boundary conditions 
on the solutions to our equations allow a steady state 
that satisfies detailed balance if and only if Local Ther- 
mal Equilibrium holds. We derive some peculiar experi- 
mentally accessible consequences along the way. 

At last, before drawing conclusions, we introduce the 
theme of a Brownian constrained particle. It was Mau- 
pertuis who first realized that a spatially dependent ki- 
netic energy could be interpreted as a deformation of the 
metric of space. Along this line, we prove that the same 
equation that describes diffusion in nonuniform temper- 
ature also describes the motion of a Brownian particle 
in m-dimensional space at uniform temperature To, that 
is constrained into a smooth curve. The equivalence be- 
tween thermodynamics and geometry suggests an inter- 
pretation of inverse temperature as a metric 



T(x) = T /£(x) 



(1) 
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where £(x) is the length of the tangent vector to the curve 
and a; is a curve parameter. In particular, the theory is 
covariant: A suitable choice of coordinates can modify 
the temperature profile, making it locally uniform. As 
a further connection between geometry and thermody- 
namics, the temperature profile can also be made glob- 
ally uniform by a coordinate transformation if and only if 
the system satisfies detailed balance. Generalization to n 
dimensions is deferred to a parallel publication [8fl , where 
the geometrical aspects are explored in greater depth. 

Before plunging into details, let us first take a glance 
at a more physical byproduct of the theory, viz. the fol- 
lowing generalization of the First Law of Thermodynam- 
ics for systems that experience a varying environmental 
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temperature (ks set to unity): 

5Q = dU-5W+(n-U/T)dT. (2) 

The 1-dimcnsional case will be given a precise formu- 
lation in the formalism of stochastic thermodynamics. 
Eq. ([2]) is a straightforward generalization for n indepen- 
dent degrees of freedom. It must be emphasized that 
both our theory and Van Kampen's imply modified First 
Laws, whose implications have yet to be fully explored, 
to our knowledge. Maybe the most remarkable fact is 
that, even keeping the internal energy fixed and perform- 
ing no mechanical work, yet heat flows when the system 
meets a temperature variation. Eq.@ is in agreement 
with the usual First Law once the last term is identified 
as a "thermal work" —SWth that needs be performed 
on the system to help it climb up against the temper- 
ature gradient, or that can be extracted when the sys- 
tem follows its natural tendency towards colder temper- 
atures. For an overdamped system, performing a transi- 
tion on the "equipartition shell" U = nT/2, the thermal 
work performed reads —SWth = ndT/2 — dU . Along an 
adiabatic transformation the total mechanical work that 
needs be performed on the system is twice the kinetic 
energy — SW = 2dU. 

2. Dynamics 

Let us make all of this more precise. In the following 
£ 4 is white Brownian noise, with (£tCf) = S(t — t'). We 
use the Stratonovich convention on the stochastic differ- 
ential, which allows to perform standard differential cal- 
culus. For sake of simplicity we omit the notation " o " 
for the stochastic differential. We contract the spatial 
derivative with d and the derivative with respect to ve- 
locities with d v . Degrees of freedom have unit inertia; 
think of a Brownian particle for definitiveness. We first 
suppose that x either belongs to a cyclic domain or to 
the real line (periodic boundary conditions), with tem- 
perature and forces smooth functions; later on we will 
comment on discontinuities and domains with boundaries 

Conventional wisdom supports the following Langevin 
equation for a particle in a bath at uniform temperature 

x t = v t , v t = -dV(x t )-jv t + v/27rCt, (3) 

whose corresponding diffusion (Kramers) equation suit- 
ably affords the Maxwell-Boltzmann distribution P* oc 
exp— T^ 1 (v 2 /2 + V) as its steady state. Generalization 
to mediums with state-dependent temperature is legiti- 
mately carried out by many authors by replacing T by 
T{x) in Eq.([3j) [3- 7], on the basis that locally, where tem- 
perature is approximately uniform, it returns the original 
equation. This is also true of any equation that con- 
tains derivatives of the temperature. Each such theory 
has its peculiar properties and leads to distinctive pre- 
dictions. In particular, the former choice leads to a dis- 
crepancy between T(x) and the notion of temperature as 



average kinetic energy, as we argue by considering the 
corresponding Kramers generator with state-dependent 
temperature 

Lvx. P = -vdP + d v [(dV + jv)P + jTd v P] . (4) 

The steady distribution in the vicinity of a given fixed 
point x turns out not to be locally Maxwellian — nor eas- 
ily computable — as one can appreciate by plugging the 
anszat R — r(x) exp — k(x)v 2 . For no choices of k(x), r(x) 
does in fact Ly.K.-R = 0. It follows that the average ki- 
netic energy of a particle at x does not coincide with 
T(x)/2. In other words, local equipartition fails. The 
steady solution only boils down to the locally Maxwellian 
form in the overdamping limit 7 — > 00, while corrections 
to local equipartition of order 7 -1 can be shown to also 
contain derivatives of the temperature [7, Supplementary 
material]. This poses an interpretational problem: In 
which sense is T(x) "temperature"? When the system is 
a binary mixture, temperature can indeed be defined as 
the average kinetic energy of the environmental compo- 
nent. Although, very often the physical nature of noise 
is not discernible, in which case temperature might be 
self-consistently defined as average kinetic energy of the 
system's degrees of freedom. In these cases Van Kam- 
pen's theory is not applicable in an obvious way. 

Alternatively, we can use the same a posteriori require- 
ment on the form of the desired Langevin equation as for 
Eq.Q, by postulating that the steady state is 

P* (v,x)oc exp - [v 2 /2T{x) + ip(x)] . (5) 

It seems reasonable in fact that locally, in a small 
region nearby x where temperature is approximately 
uniform, the velocity distribution is again Maxwellian 
and equipartition holds, irregardless of how temperature 
varies outside that small region. This assumption can in 
principle be subjected to experimental validation — al- 
though it might be complicated to control the statistics 
of particles' velocities sufficiently close to a given point. 
We expect that this assumption should hold whenever 
that small neighborhood is nevertheless big enough as to 
enable the Brownian particle to collide with particles in 
the underlying medium, all the same temperature, a suffi- 
cient number of times. In other words the mean free path 
VT/j shall be much smaller than the scale of variation 
of the temperature T/dT. 

We plug Eq.([S]) into a generic Kramers- type generator 

Lth P* = -v dP* + d v (e 2 d v P* - fP*) = 0, (6) 

where we allowed for a state-dependent diffusion coeffi- 
cient e(x), and f{x,v) is the drift. We obtain: 

-v 3 dT/2 + vT 2 dip-e 2 T + e 2 v 2 -T 2 d v f + Tfv = 0. (7) 

For this expression to vanish, the latter drift terms must 
necessarily counterbalance all others, order-by-order in 
the velocities. Hence the drift must have the form 

f(x,v) = f (x) + f 1 (x)v + f 2 (x)v 2 . (8) 
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Plugging in Eq. , after some work we obtain 

f = 3T-Tdi>, h = -e 2 /T, f 2 =d\ogVT. (9) 

The second expression is a generalized version of Ein- 
stein's relation. A convenient choice for the diffusion co- 
efficient, making us closer to Eq.©, is e 2 = 7T. It must 
be emphasized though that many other equally valuable 
models yield the same steady state. The third term 
gives a sort of Rayleight drag force, which is square in 
the velocities and depends on the temperature gradient. 
Finally we obtain the following Kramers equation for a 
generic probability density P = P(x,v,t) 

dP dP 
dt dx 



\ dP 


d 






dv 




1 dv 





and its corresponding Langevin equation 



(10) 



Vi 



= f (x t )--fv t + v 2 d log y/T(x t ) + V^T(x t ) Ct, (11) 



where /q plays the role of driving force. The left-hand 
side of the Kramers equation can be cast as (dt + L)P, 
where L is the Liouvillian 



LP 



dH(x, 7r) dP(x, v(x, 7r)) dH(x 1 tt) dP(x, v(x, n)) 



dir 



dx 



dx 



dir 



(12) 

with Hamiltonian H = T(x)ir 2 /2 + <E>, momentum tt = 
v/T(x) and $ the potential of f /T = -<9$. The effec- 
tive inertia of the particle ir/v depends on temperature, 
being smaller where temperature is higher. The steady 
probability measure in position/momentum space is the 
Gibbs measure P*dxdv = exp — H (x, ir)dx dir. 

Eqs. (|10llip are the central equations of this work. 
We will assume throughout that they locally describe 
the interaction of a Brownian particle in contact with 
a medium with varying temperature. The occurrence of 
the locally Maxwellian steady state depends on smooth 
boundary conditions. In general, different boundary con- 
ditions yield different solutions accordingly; this will be 
the subject matter of a subsequent paragraph. 

In the overdamping limit 7 — > 00, when relaxation to 
a locally Maxwellian state long precedes spatial relax- 
ation, the Kramers generator reduces to a Fokker-Planck 
generator for the spatial distribution p = J dv P. While 
Van Kampen's theory leads to the Ito generator Q , ours 
yields the Stratonovich generator 



d t p 



-0 



f p-Vrd(VTp 



= -<9J, 



(13) 



where J is the probability current. The latter diffusion 
equation corresponds to the fir st-order stochastic differ- 
ential equation x t = fo(xt) + y / 2T(xt) (t in Stratonovich 
convention. For vanishing forces, the average drift out of 
a given state x is to the hot, 



(dx t ) x = (y/2T(x t )(t)xdt 



dT(x) 



dt, 



(14) 



which marks a difference with respect to the null average 
increment in Van Kampen's theory [j| . Almost paradox- 
ically, at an equilibrium steady state particles drift aside 
even in the absence of a steady flux, a well-known phe- 
nomenon in state-dependent diffusion [9J. 
On the real line, the heat kernel 



K(x, xq, t) oc 



— , exp — — 



dx' 



(15) 



solves Eq. (|13|) with null forces and vanishing boundary 
conditions at infinity, given the point-wise initial condi- 
tion 5(x — xq) Sufficiently close to the initial state 
x — xq + Sx and at small times t ~ 5x, 



K(xq + Sx, Xq, t) 



1 (x-Xq) 

PYp ; — (16) 



to second order. Then T(xq)~ 1 gives the typical time for 
spatial relaxation near a point: The hotter, the faster. At 
larger distances and longer times there might be small 
corrections, but basically spatial relaxation is ruled by 
the value of the temperature and not by its variation, as 
far as this is not too sudden. The overdamping regime 
holds when 7 > T. The limit we posed to the valid- 
ity of the model 7 3> dT/^/T leaves enough space for 
non-overdamped regimes with 7 ~ T, as far as the tem- 
perature gradient is kept smaller than (some power of) 
the temperature, dT < T 3/2 . 

The overdamping limit is carried out in detail in Ap- 
pendix |A] Notice that, at stake with the uniform temper- 
ature case, the overdamped stochastic differential equa- 
tion cannot be obtained by just neglecting the accelera- 
tion in Eq.lfTTj). 



3. Equilibrium thermodynamic potentials 

In this and the following sections we give some phys- 
ical insights, revisiting Landauer's and Buttiker's works 
[J [ll| in the light of classical and stochastic thermody- 
namics [3, H, E3]. 

We first consider the case where the boundary condi- 
tions on Eq. pTj) allow Eq.© as the steady state. The 
first relation in Eq.® can be integrated to give 



ip(x) 



SW 



logT(x) 



(17) 



where SW — fo(x)dx is the work exerted on the system 
along displacement dx. The steady state now reads 



P* 1 
cx — exp - 



v 

2T 



5W\ 



exp-S(x,v) (18) 



where the S is an information-theoretic generalization of 
Boltzmann's entropy function of state (x, v) with respect 
to the steady state [1J| (for a motivation from Large De- 
viation Theory, see Ref. [is]). Performing the Gaussian 
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integral in Eq. ([18f with respect to velocities we obtain 
the spatial density 



p*{x) oc 



1 



: exp 



x SW 
~T~' 



(19) 



When external forces vanish, fo = 0, the steady spatial 
distribution is p cx 1/yT, reconciling to Landauer's in- 
tuitive picture. In particular, Eq. (|19p makes the current 
in the overdamped equation vanish, so that Local Ther- 
mal Equilibrium implies (global) equilibrium. If instead 
we were to impose a uniform pressure pT — const., we 
would obtain a current profile J = d log %/T from the 
cold to the hot. Analogous conclusions were reached by 
Balakrishnan and Van den Broeck in their study of the 
Jepsen gas, where one particle (the piston) is singled out 
between two one-dimensional gases of identical particles 
interacting solely by elastic collisions, and kept at differ- 
ent temperatures In that case it is the gas's density 
that settles into the Landauer state 1/y/T. In our case, 
it is tempting to figure out Brownian particles sitting on 
top of the gas as they arrange themselves according to 
the same spatial distribution as the underlying medium. 
To the notion of entropy is associated that of heat 



TdS = SQ = d (v 2 /2) -SW + dT - v 2 d log VT, 

(20) 

where we used overbraces to define the energetic and the 
thermal heat exchanges. The latter vanishes when tem- 
perature is uniform, the former accounts for the kinetic 
energy variation and for the work done by the external 
force. The first term contributing to the thermal heat can 
be interpreted in the light of Fourier's first law of heat 
conduction, and it was accounted for by Landauer. The 
second term is peculiar of our approach; it corresponds 
to the new drag force. 

The generalized First Law Eq.([2|) follows from Eq. (|2"0)) 
once the internal energy U = v 2 /2 is identified. We now 
dwell on its interpretation. The last term U/TdT is not 
unusual. In fact, from SQ = TdS = d(TS) — SdT and 
from the definition of free energy J- =14 — TS we have 



SQ = d{U - F) 



14 — F 
T 



dT, 



(21) 



wherefrom our First Law will follow, after a sensible iden- 
tification of the subsidiary thermodynamic potentials 14 
and F is done. To accomplish this we split the process 
into two subprocesses. Consider the conditional proba- 
bility P* (v, x) jp* (x) of having velocity v at a given state 
x, and define the conditional entropy Si = — \og(P* /p*) 
and the spatial entropy S2 = — logp*. The total entropy 
obviously splits into the two contributions, S = Si + S 2 - 
We define the conditional energy at a point as Ui = U. 
Hence by definition of free energy Fi = Ui — TSt we have 



F 1 (x) 



T{x) 



\ogT{x). 



(22) 



Notice that, not depending on v, at fixed x the free energy 
is a constant. Its ^-dependence suggests to take F\ also 
as definition of the spatial free energy F 2 = F\ = F, so 
that the spatial energy turns out to be 



U 2 {x) = -T{x) 



x SW 
~T~' 



With these definitions we obtain the two sub-laws 

8Q1 =dU + dT/2 - U/TdT 
SQ 2 = -SW + dT/2 



(23) 



(24a) 
(24b) 



which combined together return the above heat flux 
SQ = SQi + 6Q 2 , given U = U x + U 2 and F = 2F . 
Consider Eqs. (l2~4|) : The second is the heat balance equa- 
tion in first-order formalism, while the first provides a 
further contribution which is peculiar to second-oder the- 
ory. The average kinetic energy at a point x decreases 
exponentially to its equiparted value T(x)/2, as standard 
manipulations in stochastic calulus yield 



J t (v\ = -^[(v 2 ) x -T( X )] 



(25) 



where the average is conditioned to x at time t. In the 
overdamping limit, the relic of the internal energy dT/2 
survives m Eq. (p4b|) . while on average the first equation 
reduces to 8Q1 = dT/2; hence reductio from phase space 
to state space removes one equiparted degee of freedom 
from the definition of heat. 



4- Nonequilibrium thermodynamics 

Stochastic thermodynamics enters the game when we 
work along the solution of the Langevin equation. Plug- 
ging Eq.dlTJ) into Eq.(|2"0)): 



SQt = - [jv t - ^2 1 T{x t ) Ct \ dx t + dT(x t ). (26) 

As we abandon equilibrium thermodynamics, we shall 
consider Eq. ([26|) to be the definition of heat increment, 
on the basis that the term within square parenthesis 
corresponds to Sekimoto's celebrated definition of heat 
0, EH i with an additional thermal term which is pecu- 
liar of our approach. 

According to classical thermodynamics, equilibrium is 
characterized by SQ/T being an exact differential. In 
other words there shall exist a state function Sc — the 
Clausius entropy — such that dSc — SQ/T in all of 
its domain (condition of detailed balance, or reversibil- 
ity). A peculiarity of our definition of heat is that, when 
detailed balance holds, Clausius's entropy exists and it 
coincides with Boltzmann's. However, while the latter is 
always well-defined as minus the logarithm of the steady 
probability density, the former might not. When the tem- 
perature profile and the steady probability density are 
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smooth functions, Clausius's entropy exists if and only if 
Clausius's criterion for reversibility holds 

i ^ = -/^=0, (27) 
J(x.v) T J x T 

where the two integrals are taken along closed loops re- 
spectively in phase space and in state space. The equiva- 
lence of these two integrals is straightforward given that 
all terms involving velocities are exact differentials in 
the expression for the heat, as far as temperature is 
a smooth function. The condition in Eq. (|2"7| becomes 
more pertinent in higher dimensions, where it holds when 
/o = — TV$, with $ a smooth potential. In one dimen- 
sion a violation of Eq. (|27|) can only occur as a topological 
effect, due to a discontinuity of the external force. 

However, even a discontinuous force might lead to 
equilibrium steady states, if the probability distribution 
and/or the temperature have discontinuities themselves 
that counter-balance those of the force. This depends 
on the boundary conditions on the probability. We only 
mention, without further discussion, that a more modern 
and complete condition for equilibrium in terms of van- 
ishing macroscopic affinities [TJ EH requires to also keep 
into account discontinuities of the function TP* as an 
information-theoretic contribution to dissipation. Rather 
than a full treatment, we prefer to give an example in the 
overdamping limit to better convey these aspects and de- 
rive one experimentally accessible result. 

Consider an overdamped Brownian particle in a closed 
interval whose extremities are put in contact with heat 
reservoirs at different temperatures T c and T^. No ex- 
ternal forces are exerted. Eq. (|13p allows infinitely many 
steady solutions, for any value of the uniform steady cur- 
rent J*. Integrating d(VTp*) = -J*/VT we obtain 

fh 

J* J T-V 2 =p* c y/T B -p%y/n, (28) 

which establishes a relationship between the steady cur- 
rent and the boundary values of the probability distribu- 
tion. In Van Kampen's theory, the analogous result to 
Eq.® is 

(x h -x c )J* =p* c T c -p* h T h . (29) 

In ours, the Landauer steady state that we initially en- 
forced is the one that satisfies detailed balance, p* \JT~ C — 
V*h VTh- By means of particle reservoirs, an external 
agent can control the particle densities at the extremi- 
ties and impede relaxation, triggering directed particle 
transport. The steady probability current is given by 
the average velocity of particles by J* — J P* (x , v)v dv . 
In order to have J* ^ 0, the steady state P*(x,v) shall 
not be even in the velocities, and in particular it can't 
be Gaussian. Local Thermal Equilibrium then fails, com- 
patibly with well-known results from specific models . 
This discussion reveals that, within our theory (specified 
by the equations) , the steady state of a specific model (a 
solution to our equations) obeys Local Thermal Equilib- 
rium if and only if it is globally equilibrated. 



5. Constrained Brownian particle 

Let us now turn page. We consider a Brownian parti- 
cle living in n-dimensional Euclidean space with uniform 
temperature T . Extrinsic coordinates and velocities are 
called (z l ,u l ), with spatial derivatives <9; = d/dz l . Ein- 
stein's convention on index contraction is assumed. 

At first, the particle is assumed to obey the standard 
n-dimensional Langevin equation 

u i = f i - 1 u i + ^fFoC=-F i , (30) 

and then forced into a smooth curve with the introduc- 
tion of n — 1 independent smooth holonomic constraints 
$ a (zt) = 0. Taking the first and second time derivatives 
we obtain the following constraints on the velocity and 
the acceleration, 

s> 1 = o, e> 4 + r;; = o, (31) 

where we defined 

E? = a ( $ a , r$ = 0<0,-*°. (32) 

We assume that the matrix £" has rank n—1 and further 
define G ab = £?E$<J«, which is positive definite. 

The key step is to appeal to Gauss's principle of least 
constraint ,20} , which states that on-shell accelerations 
minimize the so-called Gauss curvature 

C(u) = 5^ - F l )(u J - F j )/2 (33) 

compatibly with the constraints. Constraints are kept 
into account by introducing Lagrange multipliers A a . Set- 
ting the first variation of Gauss's curvature with con- 
straints to zero, 

^ [C(u) + X a (E?^ + q>V)] = 0, (34) 

we obtain u 1 — F l — X a 5 i:l e" . Lagrange multipliers can be 
determined by plugging this solution into the constraint 
equation Eq. (|3"Tl) . yielding 

A„ = G a6 (r^V + Efi^). (35) 

We obtain a Langevin equation in extrinsic coordinates 

u l = -S^^G ab T b kl u k u l + (S) - G ab ^ a k ^S kl )F\ (36) 

which will now be expressed in intrinsic coordinates. We 
parametrize the curve x — > z l (x), and introduce the tan- 
gent vector t % — dz l with length 

i{x) = ^5 t3 V{x)V(x). (37) 

The intrinsic velocity Vt — x t can be expressed in terms 
of the extrinsic velocity by taking the time derivative of 
z\ = z t (xt), which yields u\ = v t t t (xt)- Hence: 

v t =£(x t )- 1 ^S ij uiui. (38) 
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Finally, we take the stochastic differential of this expres- 
sion and plug Eq. (j3"6")l in. After some standard calculus: 

v = -v 2 dlog£ + r 2 (fy - G ah Y, a k ?l b $) t l F 3 (39) 

= t~ 2 Sijt i p -jv- v 2 d\ogldt + ^2jT r 2 8 l3 f ( j , 

where we repeatedly used the constraint = 0. No- 

tice that ( t — ^ 1 8ijt' t C I 3 t is Brownian white noise. It 
is easily proven in fact that it has null average, it is 
delta correlated with unit diffusion coefficient, and it is 
normally distributed as the sum of independent normal 
stochastic variables have normal distribution. Therefore, 
identifying the force /o = £~ 2 5ijt l f 1 and, most impor- 
tantly, the temperature as square the inverse length as 
in Eq.([T]), we re-obtain Eq.(|lip. which can now be inter- 
preted geometrically as a noisy damped geodesic equa- 
tion, with the term involving Christoffcl coefficients play- 
ing the role of the Ray leight- type drag force. If the ex- 
ternal force is conservative, / 3 = —6 i: >dj&, the driving 
force /o = — Td&, whilst possibly non-conservative, sat- 
isfies detailed balance. As a curiosity, notice that in the 
expression for the square modulus of the extrinsic accel- 
eration of the particle (setting To = 1) 

a 2 = — T~ x (y 2 — v Ud log T + 4n 2 TU 2 ) , (40) 

there appears the intrinsic acceleration v and correcting 
terms, the second of which contains the squared Gauss 
curvature k 2 = 8ijdt l dP , let U — v 2 /2. 

Parametrizing the curve with the length itself, in such 
a way that the tangent vector is always normalized, 
5ijt % t 3 = 1, we obtain the standard Langevin equation 
with uniform temperature: Temperature can be made 
(locally) uniform by a (local) coordinate transformation. 
This is the tip of an iceberg: behind stands the gen- 
eral covariance of the theory under coordinate and gauge 
transformations [8[. 

6. Conclusions 

In order to sew geometry and thermodynamics, con- 
sider as an example Brownian particles floating in a fluid 
contained in a torus-shaped recipient, with minor radius 
much smaller than the major radius. Suppose to heat up 
the fluid in a sector of the torus. According to our theory, 
Brownian particles will arrange themselves in an equilib- 
rium steady state that displays Local Thermal Equilib- 
rium. Their motion is equivalent to that of particles in 
a derformed torus, with higher temperature where the 
toroidal shape contracts. Eq.([T]) says that temperature 
can be interpreted as a measure of effective length: Ran- 
dom transitions are more probable to shorter distances, 
as if due to hotter temperatures. This picture marries 
well with Smerlak's proposal of "tailoring" the diffu- 
sion of light in suitably engineered metamaterials, with 
diffusion coefficients acting as an effective curvature of 
space(time) (2lj . Several steps in our work also resound 
with the techniques used by Zulkowski and coworkers (22j 



to build geodesies in the space of thermodynamic param- 
eters out of the inverse diffusion tensor. 

Tests of our theory should be close at hand. Two re- 
cent experiments by Langon et al. @ and by Volpe et al. 
[23[ involve Brownian particles floating in inhomogeneous 
mediums. In that case all evidence suggests that the spa- 
tial inhomogeneity, due to hydrodynamic interactions of 
the fluid with walls, does not influence the steady state, 
according to a well-known theory with state-dependent 
diffusion and viscosity [24| . Similar experiments are con- 
ceivable where the medium is heated up in an nonuni- 
form way, in which case our theory predicts the spatial 
density p oc T" 1 / 2 , a Maxwellian distribution of veloc- 
ities near a given point, and prescribes how to shape a 
state space where uniform temperature Brownian motion 
occurs with the same statistical properties. Curvilinear 
Brownian motion could in principle be obtained by shap- 
ing thin recipients. This unfortunately requires to deal 
with the friction of the particle with the walls in the di- 
rection of the relevant degree of freedom. If the steady 
state profile would probably not be affected by this, all 
other statistical properties (including Fluctuation Theo- 
rems) would hardly be saved. Rather, at times when opti- 
cal traps allow spectacular control of potentials confining 
Brownian particles [25|, l26[ , it might be conceivable to de- 
sign 1-dimensional optical curves. Another experimental 
signature of our theory might regard the occurrence of an 
anomaly between overdamped and finite-inertia thermo- 
dynamics, which in the work of Celani et al. is ascribed 
to the locally non-Gaussian nature of their steady state 
[7], and might therefore be resolved in ours. 

To conclude, we used a "top-down" approach to de- 
rive an ideal theory of diffusive systems in nonuniform 
temperature environments. As far as equilibrium steady 
steady states are concerned, the theory is coherent with 
several foundational aspects in Equilibrium Statistical 
Mechanics, in particular the equivalence of statistical and 
physical entropies. Out of equilibrium, it entails a gener- 
alized version of the First Law when the system's internal 
energy and the environmental temperature differ along a 
thermodynamic transformation. The theory allows to de- 
scribe nonequilibrium particle transport due to temper- 
ature drops and particle reservoirs, prescribing testable 
relationships between boundary values of temperatures, 
densities and currents. The property of general covari- 
ance can be constructively exploited to build a geometric 
analogue. So far the environmental temperature profile 
was assumed to be given and stationary. Relaxing this 
hypothesis, it is tempting to advance a possible role of 
the theory with respect to Fourier's Law of conduction, 
whose derivation is one major open question about sys- 
tems with varying temperature 27]. "Bottom-up" mi- 
croscopic derivations are very desirable, as for example 
in the spirit of Ref. 16], with the purpose of establish- 
ing the margins of applicability of our theory and others. 
Finally, nonequilibrium stochastic thermodynamics was 
only briefly touched on. Refs. [Ill discuss Langevin- 
type formalism, but their results cannot be directly ap- 
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plied to our theory, as they don't include state-dependent 
diffusion coefficients and Christoffel terms. A suitable 
generalization of Fluctuation Theorems, Large Deviation 
properties, and the full machinery of stochastic path in- 
tegrals should and shall be applied to make precise con- 
nections between our several peculiar notions of entropy, 
heat and energy with measurable quantities. 
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Appendix A: Overdamping limit 

We derive the Fokker-Planck equation Eq. flTBI from 
the Kramers equation Eq. flT0|) . in the limit of A — » oo. 
The treatment traces out Schuss's [13, §8.2]. Solutions 
to the perturbative chain of equations are simply given, 
as they are known from previous literature. A logically 
close step-by-step derivation can be operated by means 
of expansions in terms of Hermite polynomials, as is done 
in Ref.Q. 

Let e = 7 _1 be a small parameter. We rescale time 
according to t — > e _1 s in the Kramers equation, and 
expand to order e both the generator and the probability 
density: 



tegration over the velocities, 



e 2 L 2 + eL l + L )[P + eP 1 + e 2 P 2 



0, (Al) 



where 



L = d v (v + Td v ) , 




(A2a) 


« fl ( v2dT 
L\ = —vo — o v [ 

V 2 T 


+ /0J , 


(A2b) 


L 2 = -d„. 




(A2c) 


We obtain the following expansion: 






LqPq 


= 0, 


(A3a) 


LqP\ + L\Pq 


= 0, 


(A3b) 


LqPz + L\P\ + L 2 Pq 


= 0, 


(A3c) 


The solution to the first equation is 






P {x,v,s) = q {x,s)exp 


v 2 


(A4) 


2T{x) ' 



where qo(x,s) is an undetermined function. The corre- 
sponding spatial density is found by performing the in- 



Po(x,s) = J dv P (x, v, s) = N^T(x) q (x, s). (A5) 
We also define the conditional probability 

,,2 



7r(w, s\x) 



Pq(x,v,s) _ 
p {x,s) ~ NyfTJx)^ 2T(x) 



1 v" 

exp — orTV ^ , (A6) 



which is a normal distribution in the velocity. Plugging 
P into Eq. (jA3b"l) . we obtain 

8 V (vPt + Td v P{) = v [dlog(Tq Q ) - f /T] P . (A7) 

A solution is given by 

Pi = { qi -v [dlog(Tq ) - /o/T]}P , (A8) 

where q\ = q\(x,s) is again an undetermined function. 
Substituting into Eq. (IA3c[) and integrating with respect 
to velocities gives 



dpo _ d_ 
ds ds 



P dv 



v dPidv, 



(A9) 



where it is assumed that all terms decay sufficiently fast 
at high velocities, so that 

J d v [(v + Td v ) P 2 ] dv = 0, (AlOa) 

J d v [(-v 2 dT/2T) Pi] dv = 0. (AlOb) 

Eq. (|A9| then yields 

d sPo = -Jvd{[ qi -vd log(Tq ) + vf /T} P }dv. 

(. Al1 ) 

The first integrand is odd and its contribution vanishes. 
In view of Eq. ()A5|) , we obtain 



d sPo =di [dlog(poVT) - f /T] I v 2 P dv } . (A12) 
We now perform the integration 

v 2 Pq(x, v, s) dv = po(x, s) / v 2 ir(v\x)dv 



= T(x)p Q (x,s), (A13) 
where we recognized the covariance of ir(-,x). Finally: 

^ = d {p T [d log (VTp ) - / /t] } (A14) 
= d { [VTd (VPpo) - Po/o] } (A15) 
which is Eq. (|13l) to first order in e. 
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